Positivity in the presence of initial system-environment correlation 
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The constraints imposed by the initial system-environment correlation can lead to nonpositive 
Dynamical maps. We find the conditions for positivity and complete positivity of such dynamical 
maps by using the concept of an assignment map. Any initial system-environment correlations make 
the assignment map nonpositive, while the positivity of the dynamical map depends on the interplay 
between the assignment map and the system-environment coupling. We show how this interplay 
can reveal or hide the nonpositivity of the assignment map. We discuss the role of this interplay in 
Markovian models. 
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Introduction. — The open quantum systems formalism 
is the standard tool used to understand and model the 
decoherence and thermalization of quantum systems. In 
this formalism, the total state of the system (S) and its 
environment (E), described by the density matrix p , 
evolves unitarily. However, the focus is only on the dy- 
namics of the density matrix r/^ of S by averaging the 
degrees of freedom of E. Open quantum systems are es- 
sential for physics jy , quantum information ^2] , for simu- 
lating chemistry [31 13] , and in ultrafast spectroscopy [5] . 
In many of these fields, it is customary to assume that 
at the initial time the system is uncorrelated with the 
environment. This assumption simplifies the mathemat- 
ical structure of the map. However, recently many re- 
searchers have realized that many systems of importance 
are initially correlated with the surroundings and have 
pursued investigations on systems that admit initial cor- 
relations [7] . It is well known that a system initially 
correlated with its environment may suffer from nonpos- 
itive dynamics [8 . In this article we tackle the question 
of how the initial system-environment (SE) correlation 
and the SE coupling affects the positivity of dynamics. 

The dynamical map B describes dynamics of the re- 
duced system [5HTT]. The relationship between the total 
dynamics, and the dynamics of S is shown in Fig. IT] such 
that the map is defined as the super-operator 



i5(ryo)=trB[C/po^^C/t] 



He [UA [ill] C/t] - ,1^, 



(1) 



where A is an assignment map |12H15| that captures the 
mathematical properties of the relationship between the 
reduced state and the total state. The assignment map 
captures the essence of the open quantum systems per- 
spective. It represents all the physical assumptions made 
about the total state as a function of the known reduced 
system state, containing details about the state of E and 



P'o^ 




* kavan@quantuirLlah.org 

t cesar.rodrigueziabccms.uni-bremen.de 



Figure 1. Reduced dynamics from total dynamics. The total 
state evolves unitarily pf^ — Upo^U^. The initial state of 
system (5), tjq = trsipo^), is mapped to final state rjf = 
S(?7o ) by the dynamical map B. This process may also be seen 
as ryo assigned to po by the assignment map A followed by the 
unitary transformation U (•) U'' , and finally the environment 
(E) is traced by tr^;, as in Eq. ([T|. 



SE correlation [HI [15]. The positivity of B depends 
on the interplay of the assignment map A, the details 
of the unitary evolution, and the averaging of the envi- 
ronment [THj. These three aspects cannot be isolated. 
Partial trace is a completely positive and linear opera- 
tion [17j . as is the unitary. To completely understand 
the mathematical properties of the dynamical map, the 
missing piece is to understand the role and properties of 
the assignment map. 



The assignment map was introduced as a mathematical 
mapping that takes a matrix in S* to a matrix in the SE 
space [HIS]; this is illustrated in Fig. [I] Refs. [H [TS] 
show that an assignment map is a linear, positive and 
consistent {l9j map if and only if it is of the form 
Ap [f]] = rj®T, where t is a density matrix of E (indepen- 
dent of 77) 20], i.e., it has no initial SE correlations. This 
assignment map is also completely positive, and thus the 
derived dynamical map is completely positive, indepen- 
dent of the details of the unitary. Conversely, the assign- 



ment maps for initially correlated states cannot be lin- 
ear, positive, and consistent all at the same time. Many 
researchers have examined how to relaxed the assump- 
tion of initial SE product states [TH [T31 [TBI [^Tti25] and 
have proposed physical interpretations for the nonposi- 
tivity of the dynamical map. This is important for the 
practical purpose of doing quantum process tomography 
for initially correlated SE states, see ^24, 25 . The dy- 
namical role of such correlations and nonpositive maps 
was shown to be crucial in non-Markovian dynamical 
maps [2 [5ni HZ]. Witnesses for such correlations have 
been developed ^511^. 

In this article we study the general properties of a dy- 
namical map as a function of the interplay between the 
system-environment coupling and the assignment map. 
In the real world, a system has only one particular cou- 
pling to the environment. In this paper, we focus on the 
positivity conditions when an assignment map is com- 
bined with a particular unitary evolutions and the trace. 
We begin with a brief review of assignment maps. From 
this, we find a formula to determine the positivity of the 
map that depends on the system-environment coupling 
and the assignment map. We discuss how this coupling 
can hide and reveal the nonpositivity of the dynamical 
map. We prove that it is always possible to construct a 
specific U that reveals the initial correlations by making 
B nonpositive. We also show how the coupling can hide 
the initial correlations, making the dynamics map posi- 
tive. Finally, we look at a standard class of Markovian 
dynamical models, and show how they depend funda- 
mentally on the specific couplings that hide the initial 
correlations and guarantee positivity. 

Positivity of dynamical maps. — In [14j . the relation- 
ships between SE correlations, linearity, consistency and 
positivity were summarized using assignment maps de- 
fined in terms of a set states, {Pi}, that form a matrix- 
basis for the space of S", i.e., any state of S can be written 
as a linear (but not convex) sum tjq = ^.- r. 
assignment is defined as A[riQ\ = '^^Vil 
article we will cast the assignment in a different form: 



jii J. Then the 
(g) Ti. In this 



■^Vla] = ^ akAk m ^l 



(2) 



where ak are the eigenvalues of the assignment. The 
condition of consistency is satisfied by demanding 

J^k'^kiT^ElAkVa^kl ^ ''lo- The assignment in Eq. ^ is 
equivalent to the assignments given in [lH |T5] , see Ap- 
pendix for a proof. 

The assignment takes a density matrix in the S space 
and maps it to a matrix in the SE space with correla- 
tions. For any rjQ that agrees with the SE correlations 
then A[rjQ] — pq. As a technical trick, the state of E is de- 
fined to include additional environmental degrees of free- 
dom that are not correlated with the system. Then, the 
total system-environment state becomes po = flo^^lO) (0|, 
where |0) (0| represents the degrees of the environment 
that are initially uncorrelated with the system, while Qq 
contains the correlated state. 



Lem,m,a. To generate the most general dynamics on S 
for an arbitrary assignment map, .A[?7o ] = Pp^, the total 
SE state must have the form pf^ = n^'^" (g> |0) (0|^^ 
The total space of E is split into two parts: a part that 
is correlated with S (space E^) and the remaining part 
that is uncorrelated with 5' (space Er). 

Proof. Let the action of the assignment map on 770 
yield a correlated state of SE, pQ. Now S is not corre- 
lated with anything else that it will interact with, if it 
is then we simply absorb that part into pQ. The most 
general dynamics for S then come from the most general 
dynamics of po, which is a unitary interaction with a pure 
system, see [201 133 for that proofs. We call the space of 
ilo to be SEc and the space of the pure state Er. Note 
that Pq is not a purification of ryo. It only contains the 
systems correlated to rjQ that will interact with ryo.D 

Combining Eq. (pi) with Eq. (fTl) gives 



6(r;o)-^afe(e|C/A, 77o4c/^|e) 



(3) 



The conditions for positivity for the dynamical map is 
{s\B{\r) {r\)\s) > for all {\r) , \s)} & S. That is if every 
extremal state of S is mapped to a positive operator, then 
by convexity every positive operator of S is mapped to 
a positive operator. The positivity condition in terms of 
Eq. ([3| is 

J2 ak {se\UAk\r) (r|4.[/^|se) = ^ UkWk > 0, (4) 

ek k 

where Wk = X^e I i^^W Ak\r) ^ are positive numbers. The 
positivity of B depends on the weighted sum of the eigen- 
values of A. Therefore, the values of the weights are 
important to determine the positivity of B. 

The condition for complete positivity is equivalent to 
finding the eigenvalues of B. From [S] these are found to 
be 



akz:,zr>s' {se\UAk\r') {r\Alu^\s'e) > 0, (5) 



ekrr' ss' 



where Zrs are complex numbers satisfying "^^s Ks^rs — 
1. In general this equation cannot be simplified with- 
out specific choices of A and U. Alternatively, we 
can write Eq. ([3]) as .6(770) = '^/.akBkirjo), where 
Bkirio) = tTE[UAkrioA],W] are non-trace-preserving 
completely positive super operators. Even though each 
Bk is completely positive, the corresponding ak may not 
be positive and B may or may not be completely pos- 
itive. This is because Bk are linearly independent, but 
not simultaneously diagonalizable [31] . 

What we have shown in Eqs. (|4| and (l5| is that the pos- 
itivity and complete positivity of the dynamical map are 
function of the details of the composition of the assign- 
ment map and the unitary dynamics. In the theorem be- 
low we give a mathematical construction of interactions 
U for which B is nonpositive, provided A is nonpositive. 
Then in Eq. (l6| we give a physical condition for the set of 
interactions U for which B is always completely positive. 



Theorem. For every nonpositive assignment there ex- 
ists some rj such that A[r]] = il (^ |0) (0|, where f] ^ 0. 
Then there exists a unitary transformation U, which 
leads to nonpositive dynamics for S', i.e. there exists 
\s) such that Eefe "fe {se\UAk\r) {r\Alu^se) t 0- 

Proof. — We prove this by expHcit construction of a uni- 
tary transformation violate the positivity condition given 
in Eq. Q, and therefore the condition for complete pos- 
itivity in Eq. (Is]) as well. 

First note that if the assignment is nonpositive then 
for a specific state rj the total state is not positive, and 
we have A[ri] = Q (S) \0) (0| < 0. Note ft is not pos- 
itive and therefore not a density matrix. Let us diag- 
onalize this fJ (X) |0) (0| in a separable basis [31]: ai — 
Uiil |0) (0| Ul = Y.rij \ij) {ij\ <» |0) (0| , where r^j are 
the eigenvalues of ft. 

Without loss of generality let us assume that the very 
first eigenvalue is negative rpo < 0. Although more than 
one eigenvalues can be negative, we will only need one 
negative eigenvalue. Next we have ai = roo |000) (000| -I- 
^i>o '"'oj lOjO) (OjO| + (Trcst If we take the trace with re- 
spect to E, we would get rji = f r-go + X]j ^oj] |0) (0| + 

?7rest, where ?7rest = tr£;[cri.est]- The first eigenvalue of 
rji is Too + ^ji'oj is a positive number and a^cst is a 
positive operator. Next, apply a control unitary (with 
SEc as control) that takes |OjO) to \Qji) for j > and 
leaves everything else unchanged. U2 = |00) (00| ® 1 + 

E -f r' |0j) (OjI ® vf + Y.ti' E -f 0"' m (ul ® 1, where 
v^ — EfeLo \^ ^ J) {H- The state after this transforma- 
tion is <J2 = U2CJiUl = roo 1 000) (000 1 + ro, |Ojj) (Ojj| + 
Grest- After this, apply a control unitary with E as con- 
trol t/3 = a |00) (00| + ?;/ (^ Ejf "' ® \jj) (jj| + t® 
J2i^k \J^) (J^\ ■ "^^^ state after this transformation gives 
the desired resuh. erg = C/3a2C/] = roo |000) (000| + 
'^oj \j.ij) {.i.ij\ + <^rest Taking the partial trace with re- 
spect to E we get 773 = roo |0) (0| -h J2j "^Qj b') 01 + Vrcst- 
All roj > and ?7icst is a positive operator that does not 
contain the matrix |0) (0|. And because roo < we have 
r?3<0. 

We now consider the following dynamical map from 
Eq. ([3|. We let A[ri\ = f) ^ |0) (0| and U = U3U2U1. 
This map will violate the positivity condition in Eq. (HI) 
in the main text when \s) = |0). This proves that for a 
nonpositive assignment there exists a dynamical process 
that leads to not completely positive dynamical map. D 

Pechukas [12] showed that if there are any initial cor- 
relations in SE then the assignment map is nonpositive. 
Here we have shown that the nonpositivity of this assign- 
ment map can always be revealed as nonpositive of the 
dynamics of S given an appropriate unitary transforma- 
tion. The unitary we constructed in the proof is one such 
transformation, there can be many others. 

Now that we have shown how to reveal nonpositiv- 
ity of A in the dynamics of S, we show how it can be 
hidden. For that we exploit the bipartite decomposition: 
p = rjiSiT+Xj where x is the correlations matrix [33j . Note 



that any bipartite state can be written in this form and 
trs[x] = tr£;[x] = 0. The correlation matrix has physical 
importance as it links the states of S and E. Our phys- 
ical condition and subsequent interpretation rely on this 
matrix. 

We remark that the set of unitary transformations 
{W} satisfying 



tre [W^XoVF^] = 



(6) 



lead to completely positive dynamics. This can 
be seen by noting that the action of the dynam- 
ical map is B{rjo) = tr^ [W {rjo (g) To -h Xo} W^ = 
tiE [Wtjo <»toW^ + tiE [WxoW^ . When the second 
terms is vanishing we have i3(r?o) = tr^; [Wrjo ® tq W^ , 
which is completely positive [201 ED] . 

The authors of [34] investigated the unitary transfor- 
mations that always lead to completely positive dynamics 
for any correlations; the answer turns out to be the lo- 
cal unitary transformation, U — Us Ue- This can be 
seen as a direct consequence of the Eq. ^ above since 

tvEiiUs ^ UE)xoiUs «) UeV] = tvEiUsXoUl] = 0. We 
will now see the implications of Eq. (|6| as it applies to 
models of Markovian dynamics. 

Markovian Models. — In order to highlight the signif- 
icance of Eq. m«, we will focus on its role within de- 
coherence models that rely on environmental refresh- 
ing [55H55] . A refreshing model is one where S peri- 
odically interacts with a part of E, Tn, for duration time 
T. The total state of £' is r = To ® ri (g) r2 • • • (E) t„ ® . . . . 
The SE interactions come from a unitary of the form 
Ut = exp [—itHt] where the time dependent Hamiltonian 
is^t = En^(*>T^>")K where 



e{t,T,n) 



iinT <t< {n + 1)T : 1 
for all other t : 



(7) 



and Vn is a Hamiltonian that couples 77 to t„. Fur- 
thermore, it is often assumed that each interaction Vn 
is identical to each other, except that they act on a 
different state r„. Such an unitary couples S in an 
identical fashion to different parts of E every t = riT. 
Thus, the evolution of a step of 5* is given by: jyn+i = 
tr^; [e~*-^^"?7„ (g) r„e'-^^"] = B {i]n) ■ The repeated action 
of such a map can be written as rjn+i = B" (770) . This 
is a quantum version of the Boltzmann collision model 
of the ideal gas. These models have been shown to have 
thermalization properties similar to the Markovian mas- 
ter equation for timescales much larger than T [5S1 137] . 
To understand how such a model deals with the SE 
correlations X: we will now examine the behavior of x for 
one refreshing step. Att — 0,po = r]o(^To. Thus, xo = 0. 
After coupling S and E for some time t — T, correlations 
between rji and the part of E will have developed, giving 
rise to a xi ¥= 0. However, due to the nature of the 
coupling of the refreshing model, such correlations will 
not have an impact on later steps. Note that for the next 

step, i]i will be coupled to ri, making tr£; UiXiUl = 0. 



Similarly, for each step, the correlations are discarded 
tr_E [C^nXnt^l] = 0. Eq. ^ shows how these Markovian 
models are completely positive. 

Conclusion. — We have found the conditions for positiv- 
ity for dynamical maps coming from correlated system- 
environment (SE) states. These correlations can some- 
times make the dynamical maps nonpositive, which make 
their use difficult. Thus, finding if a map is positive sim- 
plifies its use. We used linear assignment maps that can 
create SE correlations, and considered the most general 
SE couplings. Similarly, we have found the conditions 
for complete positivity of the map. 

We showed how the positivity of the map depends on 
the interplay between the assignment map and the SE 
coupling. For correlated states the assignment map can 
be non-positive, and still have a meaningful physical in- 
terpretation. The specific of the SE coupling can hide or 
reveal this non-positivity, affecting in turn the positivity 
of the dynamical map. Wc prove that if the assignment 
map has negative eigenvalues, there always exists a SE 
coupling that will reveal this negativity by making the 
dynamical map non-positive. We show how to construct 
such a coupling. 

The SE coupling can also hide the negativity of the as- 
signment map. We give an expression for the conditions 
that the SE coupling, when fulfilled, the SE correlations 
are hidden making the dynamical map completely posi- 
tive. We show how a very large class of Markovian mod- 
els, known as refreshing models and Boltzmann collision 
models, are completely positive and Markovian precisely 
because their couplings are chosen to periodically hide 
the SE correlations. 

These results highlight the dynamical role of posi- 
tive and non-positive maps in physically-motivated open 
quantum systems. This formulation explains how to use 
assignment maps to expand the dynamical map formal- 
ism to account for initial correlations and non-Markovian 
effects, expanding its utility. At the same time, these re- 
sults explains the role of system-environment correlations 
in many commonly used models. 
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APPENDIX 

The assignment presented in [T3] is of the form A[Vi] — 
Vi (E) Ti — M.i, where {Pi} form a linearly-independent 
matrix basis on the space of S, i.e. any state of S can be 
written as 77 = J2i ^i^i- 



A 



E' 



= ^r,^[P,]=^i 



(8) 



The consistency condition requires tr£;[Ri] = P^ (and 
therefore tr[Mi] = 1). Additionally, Hermiticity preserva- 



tion requires that 



. Note above, {Pj} are density 



operators but {Mi} are not necessarily positive. Here we 
show that this is the same as a map in Eq. ([2]) in the 
main text. 

Lemma 2. For any set of Linearly independent ma- 
trices {Pi}, there exists the dual set {Ai} satisfying 
tr[A, P,]=%. 



Proof. Write 



= J2j hijTj, where hij 



real 



bers and {Fj} form a Hermitian self-dual linearly in- 



dependent basis satisfying trpiFj] — 26i. 



Since 



i} form a linearly independent basis, the columns 
of matrix H — ^. hij \i) {j\ are linearly independent 
vectors, which mean H has an inverse. Let matrix 
D^ = H^^, then HD^ = I, implying that the columns 
of D are orthonormal to the columns of H. We define 
Ai = 5 X]i 'iij^j^ where dij are elements of D. D 

Lemma 3. A map in the form of Eq. (|8| is equivalent 
to the map of the form Eq. (pi in the main text. 

Proof. We write the map in Eq. (|8| as 



A[v] = Y, tr[A. v] 



(9) 



First note that by this construction Eq. ^ satisfies 
Eq. (Is]). Next, we can write the operators Mi and Ai 
in their eigenbasis; 



A[v] 



im 

E 

imn 



tr [di, 



{dimlv]/] km) {nn\ 



^ira^in yim) y^in\ ^ \^iml yi- 



(10) 

(11) 



{din\ and we 



we 



Next we define a^ = dimTin and Ak = | 
have the desired from. 

Conversely, to cast the map in the form of Eq 
have to chose a set of linearly independent matrices as 
the basis. The action of the map in Eq. (2]) in the main 
text acting on the elements of the linearly independent 

basis gives us "^i — X)fe afc^fePj^I- D 

Through out this Letter, we are use a different notation 
for assignment maps than in |14| . To aid the reader, we 
will prove that the assignment maps from [14j can always 
be written as in Eq. ([2]) in the main text. The proof is 
as follows. In [T3], the assignment map was written as 



-4W = Et^'[^j^] 



(12) 



which is clearly of the form of Eq. Q. Note that 
77, Pj and Aj are matrices in the space of S and, 
while Tj are matrices in the space of E. Note that 
tr[Aj77]Pj can be expanded using an additional index 



such that tr[Aj77]Pj 



Em lJ-m,jMm,j T] M^ ■ . AlsO, 



Tj can be expanded on its eigenbasis {\T„j)} such that 
Tj = J2n ^n.j \Tn,j) {Tn.j\i whcrc n runs up to e. Thus, 

-^W === E E m™j^".j-^^™j V mI^^ ® |T„j) (r„,,| 



This can be cast on the form of Eq. ^ in the main 
text by combining the indices k = {j, m, n} such that 

T„j). Note that A^ is 



ttfc = Mmjinj and Ak = M^j 



a rectangular matrix, mapping from S space to the SE 
space. This proves how to write Eq. (12) in the form of 
Eq. (l2| in the main text. 
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